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We report neutron scattering measurements, which reveal spin-liquid polymorphism in a ’ll’ iron 
chalcogenide superconductor, a poorly-metallic magnetic FeTe tuned towards superconductivity by 
substitution of a small amount of Tellurium with iso-electronic Sulphur. We observe liquid-like mag¬ 
netic dynamics, which is described by a competition of two phases with different local structure, 
whose relative abundance depends on temperature. One is the ferromagnetic (FM) plaquette phase 
observed in the non-superconducting FeTe, which preserves the C4 symmetry of the underlying 
square lattice and is favored at high temperatures. The other is the antiferromagnetic plaquette 
phase with broken C 4 symmetry, which emerges with doping and is predominant at low tempera¬ 
tures. These findings suggest a first-order liquid-liquid phase transition in the electronic spin system 
ofFeTei-^S.Se),. We thus discover remarkable new physics of competing spin liquid polymorphs in 
a correlated electron system approaching superconductivity. Our results facilitate an understanding 
of large swaths of recent experimental data in unconventional superconductors. 

PACS numbers: 71.27.+a 74.70.Xa 75.40.Gb 75.25.Dk 


Ideas about liquid polymorphism date as far back 
as Mendeleev’s work on water-alcohol mixtures and 
Rontgen’s proposals that many of the anomalous prop¬ 
erties of water, such as its density maximum at 4°C, can 
be understood if water is viewed as a mixture of two liq¬ 
uid ’’species”, different in local structure and density [1- 
3]. The possible existence of different liquid polymorphs 
and liquid-liquid phase transitions (LLPT) between such 
polymorphs in simple molecular fluids continues to re¬ 
ceive considerable attention, but the issue remains un¬ 
settled [4-10]. This is mainly because the competition 
between different liquid phases is governed by the en¬ 
ergy balance of inter-molecular interactions, and is thus 
expected only at low temperatures, where it outweighs 
thermal entropy. For most simple liquids this is expected 
only in a deeply supercooled region, well below the freez¬ 
ing temperature [4-10]. On the other hand, such compe¬ 
tition arises quite naturally in a system of atomic mag¬ 
netic moments in a crystal, where development of mag¬ 
netic order is hindered by competing interactions (frus¬ 
tration) [11, 12], reduced dimensionality [13, 14], or both 
[15]. In these cases, the system remains disordered even 
at temperatures much lower than the energy of magnetic 
interactions, thus realizing a spin-liquid state. Hence, 
many exotic quantum properties of liquids have been dis¬ 
covered by studying spin liquids [11-15]. 

Liquid-like states play a prominent role in theories 
of unconventional superconductivity in cupric oxides 
and iron pnictides and chalcogenides, which feature 
two-dimensional square nets of magnetic atoms with a 
strongly reduced tendency to order magnetically [16-18]. 
There is strong experimental evidence that superconduc¬ 
tivity in these materials is coupled with the antiferromag¬ 


netic correlations, which could be responsible for elec¬ 
tron pairing [18]. A simple pedestrian way to visualize 
such coupling is to imagine that an electron moving in 
a spin-liquid background, by virtue of the hybridization 
of its wave function with those of magnetic electrons, 
is ’’dressed” by a magnetic polaron, and that this mag¬ 
netic polarization cloud plays similar role to polariza¬ 
tion of atomic displacements in the phonon mechanism 
of Cooper pairing in conventional superconductors. The 
presence of delocalized electrons interacting with the sys¬ 
tem of atomic spins also affects the nature of spin liquid 
state: the emerging patterns of short-range magnetic cor¬ 
relation not only have to optimize the orbital overlap en¬ 
ergy of localized valence electrons (spin superexchange) 
but also the hybridization energy with the delocalized 
electron. 

Motivated to uncover magnetic signatures of elec¬ 
tronic pairing interactions, we carried out neutron scat¬ 
tering study of dynamical magnetic correlations in 
FeTeo.87So.13, an iron chalcogenide of ”11” superconduc¬ 
tor family, where substitution of only 13% of Te with S, 
far below any percolation phenomena, leads to the onset 
of superconductivity [19]. Through neutron scattering 
measurements, we have discovered not only the hallmark 
of spin liquid behavior, but also clear evidence of LLPT 
between two spin liquid polymorphs, one featuring FM 
plaquettes present in the parent Fei+ y Te compound [20], 
and the new one, which can be associated with the emer¬ 
gent superconductivity in FeTei_ x (S,Se)j; materials [21- 
31] (see also Supporting Information). 

The FeTei-^Sa, (x = 0.13(1), m = 2.6 g) crystal used 
in our present measurements was grown from Te-S self- 
flux using a high-temperature flux method [19] and had 
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a mosaic of < 1° full width at half maximum (FWHM). 
Elemental composition, structure and bulk properties, in¬ 
cluding the superconducting transition, of a small piece 
of the same crystal were characterized in [19]. Sam¬ 
ple was mounted on an aluminium holder attached to 
the cold head of a closed-cycle refrigerator on the HB3 
triple axis spectrometer at the High Flux Isotope Reac¬ 
tor at Oak Ridge National Laboratory, USA. Measure¬ 
ments in Figures 1 through 4 were performed in (h, k, 0) 
zone with the crystal c—axis vertical; measurements in 
the (h,0,l) zone are presented in the Supporting Infor¬ 
mation. Constant-energy scans were acquired by varying 
the wave vector transfer for fixed monitor counts in a low- 
efficiency detector between the pyrolytic graphite (PG 
(002)) monochromator and the sample and with fixed 
scattered neutron final energy E/ = 14.7 meV selected 
by PG(002) analyzer. Beam angular collimation was set 
to 48' — 60' — 80' — 240', from monochromator to detec¬ 
tor, and a polycrystalline PG filter after the sample was 
used to minimize intensity at harmonics of the desired 
wavelength. The resultant energy resolution measured 
by the FWHM of the incoherent elastic scattering sig¬ 
nal was 1.35 meV. The wave vector resolution calculated 
using standard methods (Supplementary Information) is 
shown by FWHM ellipses in Figures 1-4. The ambient 
neutron background (BG) of 9 counts per minute was 
measured with Cd-blocked analyzer entrance and sub¬ 
tracted from all data. Q -independent BG arising from 
incoherent elastic scattering in the sample subtracted in 
Figure 4 was estimated by averaging < 1% of lowest- 
intensity points. 

The experimental signatures of a liquid arise from 
well-defined microscopic local structures, fluctuating in 
time, that balance the energy of inter-molecular inter¬ 
actions and the entropy of thermal motion, which are 
at the origin of liquid state [4-10]. In a glass, another 
common amorphous state of condensed matter, which 
is solid, the local order is static and survives averaging 
over practically infinite time, giving rise to elastic scat¬ 
tering. In a liquid it is dynamic, so a well-defined local 
structure is only present as an instantaneous snapshot. 
The corresponding equal-time correlation function is ob¬ 
tained by energy-integrating the inelastic scattering sig¬ 
nal measured in neutron scattering experiment, where 
energy- and wave-vector-dependent patterns in Fourier 
space arising from time- and space-dependent real-space 
correlations are detected [32]. For a system of electronic 
spins in a crystal, it is the correlation of magnetization 
fluctuations, rather than of atomic positions, that de¬ 
termines the character of the equilibrium state. The 
hallmark of spin liquid is the presence in magnetic neu¬ 
tron scattering intensity, in a broad range of energies, of 
a well-defined diffuse pattern in the wave vector space, 
which is characteristic of a particular short-range order. 
Indeed, if such a pattern persists in the energy range 
containing the dominant portion of scattering spectral 
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FIG. 1. Dynamical magnetic correlations in FeTeo.87So.13 
measured by inelastic neutron scattering, (a)-(c) ( Q,E ) in¬ 
tensity maps along (h, 1 — h, 0) direction in the a — b plane, 
(cl), measured at 298 K, 100 K and 3.1 K, respectively. The 
corresponding (d)-(f) panels on the right show representative 
constant-E linear scans with resolution-corrected Lorentzian 
fits; the intensity is offset by adding (E-3)*30 counts to avoid 
overlaps. Dashed lines show the average peak positions, 
(g) (Q,E) intensity map, and (h) representative scans along 
(0.5, fc, 0) direction, (gl), at T = 5 K. Ellipses show the full 
width at half maximum (FWHM) of the instrument resolu¬ 
tion. 


weight, then this same pattern will dominate the single¬ 
time correlation describing the instantaneous local spin 
structure of the liquid state. 

A liquid-like, diffuse scattering pattern persisting in a 
wide energy range, appears as a broad vertical column 
of intensity in the (Q,E) coordinate plane cutting the 
4-dimensional ( Q,E) phase space along particular wave 
vector trajectory [15]. This is precisely how the mea¬ 
sured intensity is shaped in FeTeo. 87 So. 13 , Figure 1. In 
Figure 1, a-c, the ( Q,E ) slice along the (h, 1 — h) line 
of the two-dimensional (2D) Brillouin zone correspond¬ 
ing to the a — b plane of Fe square lattice is shown, as it 
evolves upon cooling from 300 K to 3 K. At h = 0.5 this 
trajectory passes through Q = the point where 

spin resonance at E ss 7 meV develops in the super¬ 
conducting FeTei-^Sea, samples [21, 22]. We observe no 
resonance even at 3 K, which is below the superconduct¬ 
ing transition temperature T c w 8 K of our material 
[19]. This is consistent with the filamentary character 
of superconductivity [19], which leaves bulk of the sam¬ 
ple probed in our neutron experiment near critical doping 
and at the verge of being superconducting. The columns 
of magnetic scattering appear at h ss 0.3 and « 0.7 in 
the units of P4/nmm (a = b ss 3.8 A) reciprocal lat- 
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FIG. 2. Electronic density and magnetization correlation in 
plaquette spin liquid models, (a)-(c) Height shows correlated 
electronic density corresponding to the isotropic atomic mag¬ 
netic form factor of Fe 2+ and enlarged by the covalence factor 
of 2, which best fits the data in Figure 3. Color represents 
the amplitude, in arbitrary units, of magnetic correlation for 
the inter-plaquette correlation length equal to two nearest- 
neighbour Fe-Fe spacings. A damped-wave correlation of four- 
iron ferromagnetic (FM) UUUU plaquettes, propagating with 
the wave vector (7r,0), gives rise to short-range checkerboard 
in FeTe, (a). For antiferromagnetic UDUD plaquettes emer¬ 
gent in FeTei_ a; (S,Se)a;, a correlation with wave vector ( 7 r, 0 ) 
gives short-range bistripe, (b), while for wave vector it 

gives collinear antiferromagnetism observed in iron pnictides, 
(c). (d)-(f) show pattern of scattering intensity correspond¬ 
ing to the model on the left, averaged over possible plaquette 
orientations to restore the macroscopic C 4 symmetry. 


tice, Figure 1, d-f. This pattern is strikingly ubiquitous 
in FeTei-^Sea; samples where superconductivity is par¬ 
tially or completely suppressed by imperfect composition, 
Fe substitution with Cu or Ni, or temperature [4-10]. 

The scattering we observe in FeTeo.s7So.13 does bear 
some similarity with parent Fei+ y Te material: it is con¬ 
sistent with similarly large Fe local moment, which is 
also obtained from the Curie-Weiss analysis of magnetic 
susceptibility [19], is stronger near Q = (n, 0), Figure 1, 
g,h, and increases in intensity upon warming, which in¬ 
dicates thermally-enhanced electronic localization. How¬ 
ever, there is also a marked distinction, which renders 
the four-spin up-up-up-up (UUUU) FM plaquette local 
structure of Figure 2a inapplicable. In Fei+^Te, where 
the overwhelming majority of magnetic spectral weight 
is also inelastic, spin-liquid-like, it can be very accurately 
described by the short-range AFM correlation of UUUU 
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FIG. 3. Neutron scattering pattern of dynamical magnetic 
correlations in the a — b plane at constant energy E = 6 meV. 
(a)-(c) data at T = 3.1 K, 100 K and 298 K, respectively, as a 
function of the in-plane wave vector ( h , k, 0 ). (d)-(e) are fits 
of the above data to the model of a mixture of two types of pla¬ 
quette spin liquids shown in Figure 2. The calculated intensity 
is orientation-averaged to restore the macroscopic C 4 symme¬ 
try observed in experiment. In addition to the background 
and the intensity scale factor, the inter-plaquette correlation 
length, £, covalent compression of the magnetic form factor, 
and the percentage fraction of the ferromagnetic UUUU pla¬ 
quettes were varied in these fits. Ellipses show the FWHM of 
the instrument resolution; dotted lines show locations of cuts 
presented in Figure 1. 


plaquettes illustrated in Figure 2, a,d [20]. The scatter¬ 
ing intensity in this model is governed by the UUUU pla¬ 
quette structure factor, which is exactly zero along the 
(h, 1 — h) line and its square-symmetric replica, giving 
rise to zero-intensity square in Figure 2, d. Hence, while 
the practical absence of scattering along these lines in 
Fei+yTe provides unique fingerprint of the FM plaquette 
liquid, this model fundamentally cannot account for the 
intensity observed in FeTeo.87So.13, Figure 1, a-c, [20]. 

The spatial structure of short-range spin-liquid cor¬ 
relation is revealed by constant-energy maps of mag¬ 
netic scattering in Figs. 3 and 4. The first important 
observation here is that all scattering features, includ¬ 
ing the characteristic ridges, which extend diagonally 
from (±0.5,0) and (0, ±0.5) positions and correspond to 
diffuse columns in Figure 1, a- e, are extremely broad. 
Resolution-corrected Lorentzian fits to the data shown in 
Figure 1, d-h, suggest correlation length of only about 
one lattice repeat, which is the distance between the di¬ 
agonal next-nearest neighbors Fe. Hence, the elementary 
local plaquette describing such correlation cannot consist 
of more than four Fe sites; the calculated structure factor 
for an ordered cluster of 8 spins already has features that 
are sharper than those seen in Figures 3 and 4. Among all 
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4-spin plaquettes, the best fit to our data is provided by 
the slanted AFM up-down-up-down (UDUD) paquette 
shown in Figure 2, b,e. In the hindsight, the presence of 
this particular local order is quite natural: its long-range 
wave-like propagation with wave vector ( 7 r, 0) describes 
bicollinear AFM partial order observed in Fei +y Te, sug¬ 
gesting that the corresponding spin liquid state, Figure 
2b, is close in energy. 

In Fei+yTe parent material, the dominant UUUU pla- 
quette correlation persists for temperatures up to at 
least T r = 300 K, and only increases in intensity upon 
warming, as a result of thermally-induced localization of 
weakly itinerant electrons (the effective spin localized on 
each Fe atom changes from S = latT = 10KtoS = 
3/2 at 300 K). That this correlation is so insensitive to 
thermal disorder suggests that formation of the UUUU 
plaquette spin liquid in Fei+^Te is governed by the hy¬ 
bridization energy notably exceeding k^T,. ?« 30 meV. 
It remained a puzzle, why partial bicollinear magnetic 
order observed in Fei +y Te below Tjv ~ 60 K, although 
it affects only minority of the spectral weight, does not 
accommodate UUUU plaquettes, thus indicating a dif¬ 
ferent liquid state. One solution to this puzzle is offered 
by the (co-)existence of two competing spin liquid poly¬ 
morphs, only one of which crystallizes into an ordered 
structure on cooling. Our present results indeed provide 
strong experimental support to this idea. A fraction of 
magnetic spectral weight in FeTeo.srSo .13 also freezes in a 
short-range glassy state below ss 50 K, corroborating that 
UDUD plaquettes are favored by energy rather than en¬ 
tropy, which favors higher-symmetry UUUU plaquettes. 

Surprisingly, the ridges characteristic of the AFM 
UDUD plaquettes in FeTeo.s7So.13 become less pro¬ 
nounced upon warming, Figure 3, b,c, while scattering 
intensity shifts to smaller Q, acquiring shape reminiscent 
of the UUUU plaquette liquid in Fei +y Te, Figure 2d. 
Similar temperature evolution is seen even more clearly in 
the pattern of (quasi-)elastic scattering, Figure 4, which 
exhibits characteristic square pattern already at 100 K. 
We therefore introduce a model where the two contribu¬ 
tions coexist, and refine the relative composition of two 
liquid components by fitting the data at different temper¬ 
atures. While this model is not very sensitive to small 
admixtures of one of the components, it clearly shows 
that that UUUU plaquettes form a minority at 3 K, but 
a majority at 300 K. That such a simple model, with only 
few parameters, provides such an accurate description of 
all experimental data (see also Supporting Information) 
is a strong evidence that it correctly captures the physics 
of the observed spin liquid state, which is a mixture of 
two components with different local structure. 

A closer inspection of scattering pattern in the UDUD 
plaquette model as a function of the inter-plaquette cor¬ 
relation length, £, leads to another surprising discovery. 
Recently, a variety of seemingly distinct, rather myste¬ 
rious and ill-understood neutron patterns, including the 
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FIG. 4. Quasielastic magnetic correlations measured at E = 0, 
showing most clearly change in the nature of dominant mag¬ 
netic plaquette with temperature, (a), (b) elastic intensity at 
4 K and 100 K, respectively. Data in (b) was measured for 
k < 0 and then symmetrized with respect to k — 0. (c), (d) 
show fits of the corresponding data on the left to the model 
of a mixture of UDUD plaquettes, (e), prevailing at low tem¬ 
perature, and the UUUU (FM) plaquettes, (f), emergent on 
warming. 


spin resonance, were observed by different groups in vari¬ 
ous FeTei-^Sea, samples, both with suppressed and well- 
developed superconductivity [21-31]. In the absence of 
even a qualitative description, these results were often 
compared with ad-hoc functions, such as a Sato-Maki 
function, leading to data parameterizations void of much 
physical meaning [22]. Remarkably, it appears that all 
the variety of the observed neutron patterns can be quali¬ 
tatively well described by the UDUD plaqutte model with 
only slightly varying inter-plaquette correlation length, 
£ ~ 1 — 2 lattice units. Hence, they all manifest very sim¬ 
ilar local correlation, revealing amazing intrinsic univer¬ 
sality uncovered by our real-space model (see Supporting 
Information). 

Thus, we find that the UDUD plaquette spin liquid 
emerges with doping in FeTei_ a: (S,Se) a; , as it begins to 
develop superconductivity. It breaks the C 4 symmetry of 
square lattice, and therefore can be related to ” nematic- 
ity” observed in many unconventional superconductors 
[18]. It is favored energetically, and, in FeTe is stabi¬ 
lized by bicollinear static antiferromagnetism. Symmet¬ 
ric UUUU plaquettes, on the other hand, are favored by 
entropy and are preferred at high temperature, as we 
observe. Our findings suggest that a liquid-liquid phase 
transition occurs in the system of electronic spins in FeTe 
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when it is doped, with Se or S, to become a superconduc¬ 
tor. The newly emerging slanted UDUD plaquettes in all 
probability reflect a new electronic hybridization pattern, 
which is likely favored by the shift of atomic energy levels 
that results from doping, and which facilitates electron 
pairing. The structure of a plaquette is thus reflective 
of particular electronic Wannier functions, which define 
the orbital content of electronic bands in the tight bind¬ 
ing view. That this change of the hybridization pattern 
belongs to the low-energy domain such that it is affected 
not only by small amounts of dopant atoms, but also 
by temperature, is a surprising discovery that calls for a 
profound revision of the tight binding paradigm. And so 
is the remarkable competition of two different electronic 
spin liquid polymorphs in a material on the threshold 
of unconventional superconductivity, which suggests new 
understanding of electronic nematicity and non-Fermi- 
liquid behavior [18]. 
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